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SELF-GENERATED  MAGNETIC  FIELDS  AND  THERMAL  FLUX  INHIBITION 


Of  criCical  concern  In  directly  driven  laser  fusion  systems  is  the 
understanding  of  electron  thermal  transport  between  the  region  where  the 
laser  light  is  absorbed  and  the  ablation  layer.  Evidence  has  accumulated 
over  the  past  several  years  which  indicates  that  the  heat  flow  rate  may  be 
strongly  inhibited.  'Die  signatures  for  strongly  inhibited  thermal  flux  are 
Increased  burn-through  times1,  reduced  Implosion  velocities2,  reduced  mass 
ablation  rates3,  density  profile  flattening1*,  significant  energy  losses  to 
fast  ions5,  and  reduced  classical  absorption6. 

Computer  hydrodynamics  models  used  to  Interpret  these  experiments  have 
typically  employed  strong  flux-limited  diffusion.  These  codes  usually  model 
the  heat  flux  as  q"1  ■  q"1  +  (fq^  )“  ^  where  ■  -KT5/,27T  is  the  classical 
thermal  flux  and  K  is  the  unmagnetized  transport  coefficient7.  q^s  - 
n  m  “l/2(kT  ) 3/2  is  the  free-streaming  flux  where  n  ,  m  ,  T  are  the  elec- 
tron  number  density,  mass  and  temperature,  respectively,  f  is  the  flux 
inhibition  factor.  Computer  modeling  of  both  planar  and  spherical  experi¬ 
ments  with  short  pulses  Ins),  high  absorbed  intensities  (>,  10 14  W/cm2)  at 
1.054  in  laser  light  has  led  to  fairly  widespread  acceptance  of  a  flux  inhi¬ 
bition  value  near  f  *  0.03.  This  approach  is  deficient  in  several  respects 
not  the  least  of  which  are  that  it  is  ad  hoc  and  it  avoids  addressing  the 
cause  of  the  flux  inhibition. 

Several  mechanisms  have  been  proposed  as  being  responsible  for  the  poor 

transport.  These  Include  ion-acoustic  turbulence8,  the  Weibel  instability9, 

large-scale  magnetic  fields  due  to  finite-sized  laser  spots  10,  small-scale 

turbulent  magnetic  fields11  and  deficiencies  in  the  modeling  of  classical 

diffusion12.  Serious  theoretical  questions  remain  concerning  the  ability  of 
Manu*cript  approved  July  11,  1984. 


Che  first  two  mechanisms  to  provide  enough  inhibition  and  it  is  unlikely 
that  magnetic  fields  generated  solely  in  the  underdense  region  or  at  the 
edge  of  the  laser  spot  would  be  effective  in  inhibiting  the  transport.  As 
yet,  none  of  these  mechanisms  have  been  fully  implemented  into  multi-dimen¬ 
sional  hydrodynamics  models  to  investigate  the  resulting  flux  inhibition  in 
a.  self-consistent  manner. 

Magnetic  fields  will  be  generated  in  laser-produced  plasmas  whenever 
the  density  and  pressure  profiles  become  noncolllnear.  13  Large  dc  magnetic 
fields  have  been  observed10  in  laser-produced  plasmas  and  have  been  asso¬ 
ciated  with  the  finite-sized  nature  of  the  laser  beam.  Modest  asymmetries 
(0(2:1))  on  the  interior  of  the  laser  beam  also  generate  large  magnetic 
fields Laser  asymmetries  generate  pressure  variations  at  the  ablation 
surface  which  provides  the  source  for  the  magnetic  fields.  As  a  result  of 
the  ablation  process,  these  fields  are  shed  from  the  ablation  layer  and 
convected  into  the  blowoff,  filling  the  overdense  region.  We  show  here  that 
these  fields  can  give  rise  to  strong  thermal  flux  inhibition.  It  should  be 
noted  that  laser  asymmetries  (0(2:1)  or  larger)  have  been  inherent  in  nearly 
all  thermal  flux  inhibition  experiments. 

A  simplified  model  equation  governing  the  development  of  the  magnetic 
field  is 

3B  *  _+  c  ♦  ♦ 

-sr  +  u  *VB  - - -  7  n  *  7  p  -  BV*u 


(1) 


where  c  is  the  speed  of  light  and  e,  ng  and  pe  are  the  electron  charge, 
number  density  and  pressure,  respectively.  The  first  term  on  the  RHS  is  the 
barodlnlc  source  term  which  generates  magnetic  fields  (and  vorticity)  when¬ 
ever  the  density  and  pressure  gradients  are  noncollinear.  Eq.  1  is  an 
approximation  to  the  complete  magnetic  field  expression7  as  obtained  from 
the  generalized  Ohm's  Law  which  also  contains  the  resistive  diffusion, 
magnetic  pressure  and  thermal-force  terms.  These  terms  have  a  very  minor 
role  in  the  simulations  presented  here.  The  resistive  diffusion  times  for 
these  fields  are  much  longer  than  the  width  of  the  laser  pulse  and  the 
magnetic  Reynold's  number  is  large  (R^  >  10),  both  of  which  indicate  that 
the  fields  are  "frozen-in".  The  magnetic  pressure  is  seldom  more  than  a  few 
percent  of  the  plasma  pressure.  At  the  ablation  layer,  where  the  fields  are 
being  generated,  the  product  of  the  electron  cyclotron  frequency  ( u>ce) 
and  the  electron-ion  collision  time  (  t^)  Is  very  small  (»ceTie  <,  10“ 2). 
This  implies  that  the  thermal-force  term15  is  not  playing  a  role  in  this 
region.  When  uceTle  >  1  in  the  blowoff  region,  the  thermal  force  term  is 
typically  an  order  of  magnitude  smaller  than  the  baroclinic  source  term.  An 
examination  of  the  complete  expression  for  the  thermal  conductivity  coeffi¬ 
cients  in  a  magnetized  plasma7  shows  that  the  ratio  of  the  coefficient  for 
transport  across  a  magnetic  field  to  the  coefficient  for  an  unmagnetized 

plasma  is  <  0.10  for  ax,  ■  1  and  Z  >  6.  This  ratio  is  -0.02  for  u>  x. 

ce  le  ce  ie 

*  4.  u  's  this  order  are  easily  attained  in  this  parameter  regime. 

We  model  the  interaction  of  moderate  to  high  Intensity  ( 10 1 3  W/cm2  <  I 
^  10 15  W/cm2),  short  pulse  (<^  1  ns),  1.054  ua  laser  light  on  thick  (25  un, 
5.5  mg/cm2)  planar  carbon  targets  with  our  hydrodynamics  model,  FAST2D.  15 
Tor  the  results  presented  here  the  nearly  sinusoidal  periodic  variation  in 


incident  intensity  is  73%  ( A1/<I>  «  0.73,  AT  ■  I  -  I  .  )  with  a  spatial 

max  nin 

wavelength  of  SO  ua  and  a  gaussian  pulse  length  of  1  ns  FWHM.  FAST2D  is  a 
fully  two-diaensional  Cartesian  code  with  a  sliding  Eulerian  grid  with  vari¬ 
able  grid  spacing.  The  ablation  layer  is  finely  resolved  with  a  grid 
spacing  of  0.10  m.  Finer  zoning  produces  no  discernable  difference  in  the 
results.  The  refined  subzoning  follows  the  ablation  front  throughout  the 
course  of  the  run.  FAST2D  solves  the  ideal  hydrodynamic  equations  using  the 
flux-corrected  transport  (FCT) 17  algorithms  with  a  two-dimensional  magnet¬ 
ized  plasma  thermal  conduction  routine  which  Includes  the  thermal  flux  per¬ 
pendicular  to  both  the  magnetic  field  and  the  temperature  gradient  (the 
Rlghl-Leduc  term)  and  the  complete  charge  dependent  magnetized  transport 
coefficient.  The  laser  energy  is  absorbed  classically  with  10%  of  the 
lrradlance  that  reaches  critical  being  absorbed  resonantly.  The  code  has 
been  extensively  documented  against  experimental  data  18  and  is  discussed  in 
some  detail  in  Ref.  16  and  references  therein. 

Figure  1  shows  the  magnetic  field  structure  at  two  different  times  for 
a  typical  1  ns  FWHM  run  (Case  E  of  Table  1).  For  this  run  the  spatially 
averaged  peak  incident  intensity  is  2.6  *  10 W/cm2  (3.6  *  10 14  W/cm2  - 
1.7  *  10 14  W/cm2)  and  the  spatially  averaged  absorbed  intensity  is 
1.2  x  10  14  W/cm2,  yielding  45%  absorption.  The  shedding  of  the  fields  from 
the  target  as  a  result  of  the  ablation  process  is  evident  at  0.4  ns.  By 
0.6  ns  this  structure  has  nearly  separated  from  the  target  and  is  being 
convected  downstream.  The  magnetic  field  structure  effectively  fills  the 
region  between  ablation  and  critical.  An  identical  calculation  except  with 
the  magnetized  transport  turned  off  gives  60%  absorption  and  peak  magnetic 
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fields  of  order  700  kG  (Case  D),  down  from  3MG  in  the  previous  case 
numerical  results  are  summarized  in  Table  I. 

TABLE  I 

Data  Comparison  at  Peak  of  Laser  Pulse 


47.0 

143.0 

0.3 


1.6 

1. 

2. 

0. 
0.94 
24.0 
16.0 
156.0 


9.0 

25.0 

25.0 

7.2 

16.0 

12.0 

9.9 

22.8 

12.8 

2.3 

3.5 

1.5 

1.7 

5.0 

2.8 

6.0 

3.0 

17.0 

144.0 

236.0 

117.0 

127.0 

190.0 

400.0 

0.73 

0.67 

3.0 

M(U)  are  cases  with  the  magnetized  (unmagnetized)  transport. 

<I>2  (<I>A)  is  the  spatially  averaged  Incident  (absorbed)  laser  intensity 
(1013W/cm2);  P(Mb)  is  ablation  pressure;  m  (10s  gm/cm^)  is  mass  ablation 
rate;  v(106cm/s)  is  target  velocity;  Ap/p  is  variation  in  ablation  pressure 
across  target;  Dac  ( un)  is  distance  between  ablation  and  critical;  L(  nn)  is 
density  scalelength  at  1/4  critical;  B(MG)  is  peak  magnetic  field  strength. 


The  magnetic  fields  Influence  the  evolution  of  the  laser  plasma  Inter¬ 
action  In  several  ways.  The  fields  are  generated  at  the  ablation  surface 
early  In  the  pulse  near  where  the  gradient  In  the  laser  Intensity  Is  the 
largest;  approximately  midway  between  the  Intensity  maximum  and  minimum. 

The  fields  are  shed  from  the  ablation  layer  and  reduce  the  longitudinal 
transport.  The  Righl-Leduc  term,  which  controls  the  flux  perpendicular  to 
both  B  and  71,  plays  a  dual  role  here.  This  term  produces  a  transverse  heat 
flow  from  the  cooler  to  hotter  regions  in  the  underdense  plasma.  This  Is 
the  source  for  the  thermal  instability  discussed  in  the  literature. 19  For¬ 
tunately  this  Instability  is  self-limiting.  An  increase  in  the  transverse 
temperature  gradient  in  the  underdense  magnetized  plasma  also  increases  the 
rate  of  heat  flux  out  of  the  system  back  towards  the  laser.  Also,  an  in¬ 
crease  in  temperature  in  the  region  where  the  laser  is  hottest  reduces  the 
amount  of  classical  absorption  there  and  effectively  turns  off  the  driving 
term  for  this  instability.  The  temperature  throughout  the  underdense  region 
is  increased  as  a  result  of  the  inhibited  transport  thus  reducing  the  amount 
of  classical  absorption. 

Figure  2  compares  the  density  (p),  pressure  (p)  and  temperature  (T) 
profiles  which  have  been  spatially  averaged  across  the  system  for  magnetized 
transport  and  unmagnetized  transport  for  approximately  the  same  absorbed 
intensity.  Relative  to  the  unmagnetized  case,  the  profiles  for  the  magne¬ 
tized  case  show  a  reduced  ablation  pressure,  steeper  temperature  gradient  and 
higher  temperature  in  the  blowoff  and  flatter  density  profile. 


The  mass  ablation  rate  Is  very  sensitive  to  the  thermal  transport. 

Fig.  3  compares  the  mass-ablation  rates  between  the  inhibited  and  uninhi¬ 
bited  cases  as  a  function  of  average  absorbed  intensity.  For  the  uninhib¬ 
ited  transport  case  m  which  agrees  well  with  experiment  and  theory 

in  the  low  intensity  regime20  (I  <  2  -  3  *  10  13  W/cm2)  where  the  transport 
is  not  expected  to  be  strongly  inhibited.  For  the  inhibited  transport  case 
the  mass  ablation  rate  is  a  strong  function  of  the  absorbed  intensity  and 
shows  two  fairly  distinct  scaling  regimes:  m  a  1*^  for  I  <  10 14  W/cm2,  and 
m  a  I*  for  I  >  10 14  W/co2.  Also  shown  in  Fig.  3  are  the  mass-ablation 
rates  from  a  1 -dimensional  simulation  (FAST1D)21  with  an  imposed  flux 
inhibition  factor  of  0.03.  The  agreement  is  quite  good. 

In  summary,  we  have  shown  that  the  observed  thermal  flux  inhibition  can 
stem  directly  from  the  strong  magnetic  fields  generated  at  the  ablation 
layer  as  a  result  of  modest  laser  asymmetries.  With  one  notable  excep¬ 
tion22,  where  little,  if  any,  flux  inhibition  was  observed,  laser  asymme¬ 
tries  have  been  inherent  in  nearly  all  flux  inhibition  experiments.  The 
multlmegagauss  magnetic  fields  are  shed  from  the  ablation  layer  and  fill  the 
overdense  region  strongly  influencing  the  thermal  transport.  The  self-con¬ 
sistent  treatment  of  thermal  transport  in  this  environment  shows  strong 
thermal  flux  inhibition  which  is  manifested  by  the  following  observations: 
(1)  reduced  ablation  pressures,  (2)  reduced  implosion  velocities, 

(3)  reduced  mass  ablation  rates,  (4)  density  profile  flattening  and 
(5)  reduced  classical  absorption;  all  of  which  have  been  experimentally 
observed.  The  spatially  averaged  mass  ablation  rates  obtained  from  the 
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Figure  3 

Mass  ablation  rate  (gm-cm2/s)  as  a  function  of  the  spatially  averaged 
absorbed  intensity  (W/cm2)  for  the  unmagnetized  (•)  and  magnetized  (x) 
transport  and  the  one-dimensional  results  with  f  ■  0.03  (0).  A  is  the 
result  for  1/2  ns  FWHM  pulse  and  |  is  the  result  for  a  1/4  ns  FWHM 
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